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I Introduction 



a) General considerations: 

The formulation of the Ashtekar variables [l|] has led to a considerable body of work 
applying them to problems of quantum gravity. These variables are the result of a 
complex canonical transformation on a set of 3+1 Hamiltonian variables of general 
relativity related to the ADM variables. One point about canonical transformations 
that is perhaps slighted is that they sometimes allow one to find particular quantum 
solutions. The usual mapping of quantum solutions from one set of canonical variables 
to another is generated by multiplying a solution in one set of variables, by e lG , 
where G is the generator of the particular canonical transformation to give ip = e tG ^> . 
Notice that if one manages by some technique to find a particular solution \l/ in 
one set of variables and the generating function is known, it is possible to find a 
more complicated solution i[) in terms of the old variables. This technique has been 
used with the Ashtekar variables to find a few exact solutions in quantum cosmology 

ill- 

There are a number of ways in which this concept can be extended beyond these 
results. One is to attempt to promote quantum cosmology solutions to full quantum 
gravity solutions as was done with the Chern-Simons solution in terms of the Ashtekar 
variables^]. Another route is to study the concept of canonical transformations on 
the usual 3+1 variables and attempt to use the new variables to generate new exact 
solutions to quantum gravity. 

The plan of this article is to investigate several canonical transformations related 
to the transformation that gives the Ashtekar variables and use them is the sense 
mentioned above to generate new exact solutions. At this point all the solutions we 
have been able to find are minisuperspace (quantum cosmology) solutions, specifically 
solutions for diagonal Bianchi type IX cosmological models (Mixmaster models). We 
will begin with a brief discussion of canonical transformations in Hamiltonian formu- 
lations of general relativity and their relation to quantum gravity. We will then use 
these concepts to write down a series of equations that can be solved to give families 
of exact solutions. Next we will discuss the form and meaning of these solutions. The 
last section of the article will be devoted to the significance of our solutions in the 
context of quantum gravity and suggestions for using canonical transformations to 
find new solutions in quantum gravity. 

b) Canonical transformations and Ashtekar Variables: 

Both the ADM Hamiltonian variables for the gravitational field and the Ashtekar 
variables are based on a 3+1 decomposition of the four-metric of space time. The 
ADM variables consist of the three-metric components g a b on t = constant slices, 
and their conjugate momenta ir ab constructed from the extrinsic curvature of these 
surfaces, K ab . The Ashtekar variables result from a "Bargmannization" ]5| of these 
variables that is based on a complex canonical transformation similar to p — > f(q)+ip, 
q — > q in an ordinary one-dimensional classical system. It is useful to look at a similar 
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transformation in terms of the ADM variables. If we write the ADM action as 

/ = / {K ab 9ab ~ N (^G abcd 7r ab n cd -^(( 3 R- 2A) ) j + N a -K a \ b }d\ (1.1) 

where, as usual, 3 R is the Ricci scalar on t = const, surfaces, | is a covariant derivative 
on these surfaces, A is the cosmological constant, and G a bcd is the DeWitt metric, 
Gabcd = gacQbd — (1 /2)g a bg c d- The parameter ( is used to control the signature of 
space-time; it is chosen to be +1 for Lorentzian signatures and -1 for Euclidean 
signatures. It is possible to attempt the complex canonical transformation 

9ab -* 9ab, (1.2a) 

n ab -> r ab , r ab = f ab {g cd ) + in ab . (1.26) 

Here we assume that f ab could depend explicitly on g ab)C as well as g ab . The trans- 
formation is canonical (ignoring topological complications in the space of metrics) if 
and only if 

5f ab 5f cd , , 

• =0. (1.3) 



89cd Sg. 
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The action now becomes 



1 = I { (1) Tah ' 9ab + ^ fab ^ ab) - N (^l- G ^r ab r cd + 2G abcd f ab r cd - 

-G abcd f ab f cd ] - ^(C 3 R - 2A)) - 2iN a (r ab lb - f ab \ b )}d A x. (1.4) 

Notice that the term if ab g a b is a total time derivative because (1.3) implies that f ab 
is the functional derivative of a functional S, i.e. f ab = 5S/6g a b- It is also possible to 
remove the three-curvature term from the Hamiltonian constraint if we take 

G abcd ^^ + g( 3 R = 0, (1.5) 
og a b og a b 

which is the Einstein-Hamilton- Jacobi equation |J for S with the sign of g 3 R reversed. 

There are a number of observations that we can make about the action (1.4). First 
of all, if S is taken to satisfy (1.5), the Hamiltonian constraint, 

H ± = = -G abcd r ab r cd + 2G abcd -^r cd + 2gA, (1.6) 

og a b 

is an algebraic function of r ab with at most second-order terms. An obvious solution 
to 7i± = (when A = 0) is r ab = 0, which is the complex equivalent of the usual 
Einstein-Hamilton- Jacobi formulation where from (1.2b) 7t ab = i5S/5g ab , S obeys 
(1.5), and the space constraint reduces to (SS/5g ab )\ b = 0. 

Notice that the complex version of the canonical transformation is not necessary, 
one can define n ab — > r ab = 6S/5g ab + 7r ab , and the action becomes 



r ab g ab - N(^-[G abcd r ab r cd - 2G abcd —r cd + 

V9 dab 
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+G abcd ^-^-} - ^(C 3 R - 2A)) + 2N a [r ab {b + (^-) }}d\ (1.7) 
og a bog C d 1 \og a b)\ b 

and if 5 obeys the ordinary Einstein-Hamiltonian-Jacobi equation 

G abcd ^^-g( 3 R = 0, (1.8) 

r afc = is a solution to the Hamiltonian constraint (for A = 0) and (SS/5g ab )\ b = is 

again the content of the space constraint. Of course it is not necessary to assume that 
T ab ig 

equal to zero, or that S be a solution to the Einstein-Hamilton- Jacobi equation. 
If S is assumed to be any function of g a b (and its derivatives), then, for example, in 
(1.7) the Hamiltonian constraint becomes 

H ± = -^ F [G abcd T ab r cd - 2G abcd ^r cd ] - VsC (^' - 2A), (1.9) 

Vg bgab 

where 3 R' is a new "scalar curvature" defined by 

3 -R = 3 R — —Gabcd-, — t — • (1-10) 
(g dg ab dg cd 

In view of the fact that there are a number of possible linear combinations of 
n ab and SS/5g a b that are, in principle, acceptable, we would like to study general 
transformations of the form 

n ab _^ T ab = SS /5g ab + (3n ab , (1.10a) 

g ab ^g a \ (i-iofe) 

where S will not be assumed to be a solution of the resulting Hamilton- Jacobi equation 
and r ab will not necessarily be taken to be zero. Which of the approaches outlined 
above one chooses depends on the system one is studying and the goal one is trying 
to achieve. For the classical theory it might seem to be less desirable to use an S 
that does not remove the curvature term from the Hamiltonian constraint, but if one 
is using a complex canonical transformation, the quantum theory can be made more 
difficult by the necessity of imposing a reality condition on quantum states, and one 
can trade the existence of 3 R ^ for explicitly real quantum variables 0. 

Notice that even for complex canonical transformation we mentioned above there 
are, in principle, as many such transformations as there are solutions S to (1.5). The 
main difficulty in finding functions, S, of g ab and g a b,c is that f ab = 5S/5g a b is a 
two-index object, and it is difficult to construct such an object that satisfies (1.5) 
solely from g ab and g abjC (or the Christoffel symbols T bc ). However, if one considers an 
orthonormal basis of one-forms on t = const, surfaces, a 1 = e a (x c )dx a , (ds 2 = o % a l \ 
then the connection coefficients T l - k have the natural symmetry = and it is 
possible to construct the spin coefficients = {l/2)e' llk T(kj. Using these, Ashtekar 
was able to find an elegant solution for the equivalent of S. 

Before we continue, we introduce some notation that will be used throughout the 
paper. SO (3) indices will be denoted by lower case latin letters from the middle of the 
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alphabet, i, k, £ • • • — 1, 2, 3 (we reserve letters from the beginning of the alphabet 
for tangent space indices). We will use additional indices /, J, K, ■ ■ ■ — 1, 2, 3 as 
labels for certain geometrical objects. The 3-dimensional Levi-Civita tensor density 
and its inverse are denoted as fj abc and r] a bc and the internal Levi-Civita tensors will be 
denoted as €uk, Sijk- The basic fields in the ADM formalism with an 5*0(3) internal 
symmetry are a densitized triad Ef and an object K l a closely related to the extrinsic 
curvature. We introduce e\ as the inverse of e l a (the coefficients of a' 1 ) satisfying 
e i e l = Ki e i e i = $1- The determinant of e l a is defined as, 

detei = e = I^% fc e^. (1.11) 

Finally, the S0(3) connection T l a (where Y l a = r*-e£, defined above) compatible 
with e\ is 

ri = ~Y~(<4 - 2eiel)fj bcd d c e dj . (1.12) 

The Ashtekar variables are a densitized basis Ef = 2ee" and the equivalent of r ab , 
A\ = Y\ + iK' l a . We will not go into the details of of the ADM action written in terms 
of these variables (although we will give the constraints later) but just point out that 
Ef and A a { are new canonical variables derived from Ef and K l a by means of a complex 
canonical transformation generated by the equivalent of S, S — 2i J EfT l a d 3 x. Notice 
that we have an explicit form for S, whereas in the previous formulation embodied in 
(1.2) we would have had to find a solution to the Einstein-Hamilton- Jacobi equation 
in order to give an explicit expression for S. Introducing Ef as a basic variable 
introduces a new symmetry, the freedom to perform SO (3) rotations in the ijk-mdices 
without changing g a b. This symmetry is mirrored in a new constraint which needs to 
be added to the usual diffeomorphism and Hamiltonian constraints. The constraint 
structure in the new variables is 

V a ^ a = 0, (1.13a) 

K b Ef = 0, (1.136) 

e ijk EfE b F abk + 2(detEf)A = 0, (1.13c) 

where (1.13a) is the basis-rotation or Gauss-law constraint, (1.13b) is the diffeomor- 
phism constraint, and (1.13c) is the Hamiltonian constraint (modulo a term propor- 
tional to the Gauss law constraint). The quantity 

Fi lb = 2d [a ai ] +e%AiA k b , (1.14) 

is the curvature of the connection A l a , and V a Aj = d a \i + e l j k A a j\k is the SO (3) 
covariant derivative acting on internal indices. 

Of course, the definition A\ = T l a + iK\ is equivalent to r ab = SS/Sg ab — iir ab , 
and as we mentioned above, the possibility of other linear combinations of momenta 
and functions of the metric (here triad) still hold. It should be possible to find 
a generator for a real canonical transformation of the form A\ = — f l a + K l a that 
removes the potential term from the Hamlitonian constraint just as the transformation 
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T ab = —5S/5g ab + 7c ab does if the generator S satisfies the equivalent of the normal 
Einstein-Hamilton- Jacobi equation. Unfortunately, no such generator is known at 
this time. It is also possible to construct general transformations of the type given by 
(1.10) which do not remove the curvature terms from (1.3). In fact, one of us (F. B.) 
has studied canonical transformations of the type Ef = E" 1 , A % a = Y l a + f3K l a where j3 
was taken to be -1 and the Hamiltonian constraint became 0] 

e^ k E?E'(F abk - 2R abk ) - (detEf)A = 0, (1.15) 

where R l ab is the curvature of the connection Y l a) 

Kb = 2d [a r b] + e) k T j a T k b . (1.16) 

As mentioned above, real transformations have the advantage (at the possible cost 
of more complicated equations) of not requiring reality conditions in the quantum 
formulation. 



c) Minisuperspace Models 

In order to give concrete examples of possible Hamiltonian formulations that we 
discussed above, we would like to apply them to minisuperspace models where, as 
we will show, they lead to a number of exact solutions in the minisuperspace sector, 
some of which were known and some of which are new. 

The minisuperspace examples we will use are the diagonal Class A Bianchi cos- 
mological models where the metric has the form 

ds 2 = -dt 2 + 5u (t)crV J , (1.17) 

where gjj is a diagonal matrix and the a 1 are invariant one-forms that satisfy da 1 = 
(l/2)C I JK a J Aa k where the C 1 JK are structure constants of the form Cj K = ra IL e LJK , 
where m IJ is a matrix of constants ||. We will be most interested in the Bianchi IX 
case where m IJ = 5 IJ . 

For Class A Bianchi models we can write the ADM action in terms of gu{t) and 
the basis components of the ADM momentum as 

I = J W J 9u - N (^[G ijkl it ij tt kl } - ^(C 3 R - 2A)) jdta 1 A a 2 A a 3 , (1.18) 

since for all vacuum diagonal Class A models the diffeomorphism constraint is iden- 
tically zero. Here 3 R is an algebraic function of gu{t) and the structure constants 
C J JK (or m IJ ), and g = det(g u ). 

It is now possible to make the same canonical transformation t ij = f IJ + m IJ 
as in (1.2), but the advantage in the minisuperspace is that we can realize f IJ as 
dS/dgu, replacing the functional derivative by a partial derivative. The Einstein- 
Hamilton- Jacobi equation (1.5) becomes a partial differential equation, 

Gijkl-r — -r — + g( 3 R = 0. (1.19) 
dgu dg KL 
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This equation has a number of particular solutions for Bianchi IX models, some of 
which have been given elsewhere [|3| || and several more which will be discussed below. 
Notice again that this equation contains the entire classical problem (both Lorentzian 
and Euclidean), so, in principle, it has a rich solution space and consequently a large 
family of canonical transformations of the form r IJ = (dS/dgu) + iir IJ . 
As before, t ij = is a solution to the Hamiltonian constraint (A = 0) 

H ± = = -G ijkl t ij t kl + 2G ijkl —t kl + 2gK, (1.20) 

dgu 

but this is not the only possible solution, and as before, linear canonical transfor- 
mations of the form t ij = (dS/dgu) + f3n IJ are possible, where S may or may 
not be chosen to annihilate the 3 R term, we will discuss these possibilities in the 
minsuperspace context below. 

While we will refer to the connection between the formulation in terms of the S 
solutions given above and the Ashtekar variables and similar variables, we will gen- 
erally work in terms of the variables themselves in order to write the minisuperspace 
quantum equations. Writing the metric in the form of (1.7), the one-forms a 1 are 
a 1 = a I a (x c )dx a , and the orthonormal one- forms e l a are e\{t)a I a . The variables given 
above become: 

e\ = e\(t)ai(x c ), (1.21a) 

K = 4(*KV), (i-2i6) 
K = *4(*KV), (i.2ic) 

El = £/(t)det(aK(x c ), (1.21c?) 

where all the x-dependence is contained in a^(x c ). The basic variables for our pre- 
sentation will be functions of t (and m IJ ) only that can be written in the form 

F] k = {a\m IL e LJK + 4 s a>|,)e/ e f , (1.22) 
1 [5\e j je jK m KJ - 2e\e lK m IK ], (1.23) 



t 2(dete}) 
and the curvature 

e^ k E<tE)R ahk = -2(det^) 2 e^ 4j ei,e Lj (m iX m 7J - 2m LJ m IK )], (1.24) 

which, since it is a density, retains the x-dependent determinant (detcr^). 

Diagonal Class A Bianchi models, as we will see below, satisfy identically the 
Gauss-law and diffeomorphism constraints, so the only constraint that survives and 
provides the quantum operator we will need to determine the minisuperspace wave 
function of our models is the Hamiltonian constraint. 

d) The Quantum Problem 

As we mentioned above, we would like to use quantum minisuperspaces as models 
of quantum gravity in which it is possible to find exact particular solutions that can 
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be used to compare quantization in the different sets of canonical variables discussed 
above. We will use the Dirac scheme of quantization where we will apply the Hamil- 
tonian constraint to a state function \l/ and obtain a form of the Wheeler-DeWitt 
equation for the function In order to do this we will realize some of the operators 
in the quantum system as derivative operators on functions of the others. In the ADM 
formulation, for example, the metric gu and its conjugate momentum tt ij become 
operators and one can make ^ a function of gu and realize ir IJ as —id/dgjj. It is 
also possible to choose the "momentum representation" in which \1/ is a function of 
7r /J and gu is realized as —id/dn IJ . It is also possible to have ^f(gu) and realize 
t ij as d/dgu as in the Bargmann-Segal formulation || or choose the "connection 
representation" in which \1/ = \l/(r /J ) and gu becomes —d/dr IJ . Since a} in the 
Ashtekar representation is the conjugate of E( , we have the same possiblities, that of 
E- — > E-, a\ —> djdE\ or E\ — > d/da}, a) — > a). We will investigate a number of 
them. 

Notice that for any solution to (1.19) the operator version of (1.20) in the tjj — > 
d/dgu representation (with A = and all derivatives standing to the right) has the 
form 

d 2 OS d 

T~L± — —Gu KL — \-2Guklt; — t; , (1.25) 

ogijdg KL dgu dg KL 

which has as a solution to 7ijJ& = 0, \I/ = — const. By the usual transformation 
of variables we can construct a solution of the usual Wheeler-DeWitt equation of the 
form ip = e ±s ^o. This formal solution becomes a true solution if we have an explicit 
solution for S. Notice that for the usual Hamilton- Jacobi formulation a solution 
ip = e ±lS ^o is possible if a solution S can be found for the usual Einstein- Hamilton- 
Jacobi equation. It is also possible to study the quantum problem for any of the linear 
canonical transformations of the form gu — ► gu, tt ij — > f IJ + /3tt ij . Of course, the 
equivalent Ashtekar-like transformations also lead to quantum equations and there 
exist similar maps among the quantum solutions. Our plan is to investigate a number 
of these possibilties, present several exact solutions, and use them as model examples 
of possible quantum solutions in the full theory of gravity and discuss the relation 
between them and such solutions. 



II Solutions to the Hamilton— Jacobi Equation 

The purpose of this section is to discuss new solutions to the Hamilton- Jacobi equa- 
tions as an intermediate step to finding solutions to the Wheeler-DeWitt equation 
for Bianchi IX models. We do this both in the real ADM and the real Ashtekar for- 
mulations. Whereas in the ADM case the theory is explictly real for both Euclidean 
and Lorentzian signatures (that we describe collectively by using the parameter ( 
introduced above) the usual way to treat Lorentzian signatures with Ashtekar vari- 
ables is by working with complex fields and imposing some "reality" conditions that 
can be used, for example, as a tool for fixing the scalar product. A less conventional 
attitude is to use explicitly real "Ashtekar-like" variables (in the sense that they keep 
their geometrical meaning)^. The Lorentzian theory is recovered by modifying the 
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Hamiltonian constraint through the introduction of a potential term. In this paper 
we choose to concentrate on this second, and more novel, approach, motivated by the 
desire to know whether this new formulation provides us with a useful alternative to 
the use of reality conditions. To this end it proves to be convenient to compare the 
results obtained with ones corresponding to the ADM case, so we begin by studying 
the geometrodynamical formulation. In order to facilitate the comparison between 
the ADM and Ashtekar formalisms we slightly modify the ADM constraints by in- 
troducing an internal SO (3) symmetry. For the specific example of the Bianchi IX 
model we get the constraints 



e ljk k j jE kI = 

e'j^Ei - 2e l]k e\e Nk 5 IN kl = 

(detaf {2C [(Tre 2 ) 2 - 2(Tre 4 )] - 2A(detE I i ) + 2k 1 } k^} £?/£?/} = 0, 



(2.1) 



where Tre 2 = e^ejO -7 " '5ij and Tre 4 = e l I e'je k K e l L 5 I J 5 K L 5 j k 5u and ( controls the space- 
time signature (£ = ±1 for Lorentzian and Euclidean signatures respectively). In 
this paper we concentrate on the study of Mixmaster models for which k} and E{ are 
taken to be diagonal, 





' V 





" 




' M 





" 


k\ = 





V 





, El = 





N 













A 










L 



(2.2) 

the variables //, v, A, M, N, L introduced above are canonically conjugate pairs, i.e. 

(2.3) 



{M,n} = l, 
{N,u} = l, 
{T,A} = 1, 



and the remaining Poisson brackets are zero. The previous expressions, together with 
(1.21d) allow us to write 



1 



V2 



(LMN) 1/2 



l/M 
l/N 
1/L 



(2.4) 



Introducing (2.2) and (2.4) in (2.1) we find that the first two constraint equations are 
identically satisfied and the scalar constraint is given by 



2(fxuMN + fxXML + v\NL) + ((M 2 + N 2 + L 2 )- 

-^ 2 ( M ^ + ^ + ^f)-2AMNL = 0. 



(2.5) 



The usual Ashtekar constraints for type A Bianchi models are (m IJ = 5 IJ gives 
Bianchi IX) 

e l]k a{E u = 0, (2.6a) 



{m IL e LJK a\ + e ijk ajj a Kk )E^ 



0, 



(2.66) 
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e ijk \a\m 1L e LJK + e l lm a j ja k K \ EfE^ = 0. (2.6c) 

With this form of the constraints real variables describe Euclidean signatures and 
complex variables, with the addition of reality conditions, Lorentzian signatures. For 
Bianchi IX diagonal models we write 



(2.7) 



where the pairs (a, A), ((3,B), (j,G) are canonically conjugate Q. If we use the real 
formulation in terms of Ashtekar-like variables introduced in Ref. pj and follow the 
same steps as before we find that the Hamiltonian constraint has a potential term. 
For the Bianchi IX model the constraint is given by 





a 





" 




' A 








4 = 





p 





, El = 





B 







_ 





7 . 










G 



'-ijk 



1L 



£ljk + e i m aja K 



+4(e i I e iJ ){e :i K e Lj ] 



2A(det£); 
0, 



{21 



m LK m IJ — 2m LJ m IK 

with m IJ = 5 IJ . If we introduce (2.7), the Gauss law and vector constraints are 
identically satisfied (as before) and the scalar constraint (2.6c) becomes 

2(aP + j)AB + 2(a 7 + 0)AG + 2(/5 7 + a)BG + 2A ABG = 0, (2.9) 

while the constraint (2.8) for the real Ashtekar-like variables becomes 



2(ap + i)AB + 2(a 7 + (3) AG + 2(/3 7 + a)BG+ 
+2(A 2 + B* + C 2 ) - + ^ + ^W-) ~ 2A ABG 



0. 



(2.10) 



In conclusion, we see that to quantize these models we only need to consider the 
Hamiltonian constraints (2.5), (2.9), and (2.10). In all the cases we can use either a 
"position" representation or a momentum representation. Since both (2.5) and (2.10) 
are non-polynomial in some of the variables, we multiply them by appropriate factors 
in order to avoid the appearance of derivatives in the denominators in some of the 
quantizations of the model. The unpleasant consequence of this is that the differential 
equations that the wave functions will have to satisfy are of very high order. 

Starting from (2.5), if we quantize by realizing the operators /}, z>, A, M, N, L as 



p>,i>,\—> 
M, N, L 



fx, v, A, 



id \, 



(2.11) 



(ft = 1) we have [/}, M] = -i, [z>, N] = -i, and [A, L] = -i. The Wheeler-DeWitt 
equation now becomes (choosing the operator ordering corresponding to writing all 
the derivatives to the right and multiplying the constraint by N 2 M 2 L 2 ) 



{2 [iavd,d u + /iAc^A + v\d v d x ] dldldl + CdldldHdl + % + d\)- 
-^(dffi + dpi + &t&i) - 2zA^f} ^ = 0. 



(2.12) 



1 For the Euclidean theory this means just that {a, A} = {f3, B} — {7, G} — 1 with the remaining 
Poisson brackets equal to zero. If we consider complex variables the previous Poisson brackets pick 
a purely imaginary factor and become {a, ^4]} = B} = {7, G} = i. 
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If instead, we quantize using 



A — > -id M , -id N , ~id L , 
\I. X. I. -> M, AT, L, 



(2.13) 



we get 



{-2(MNd M d N + MLd M d L + NLd N d L ) + ({M 2 + N 2 + L 2 )- 

-K + ^ + ^T) - 2A MiVL} * = 0. (2 - 14) 

We consider now the Hamiltonian constraint (2.9). If we quantize according to 

Pi 1 * ®-i Pi 1 ' i ^2 \^>) 

A,B,G-^ —id a , —idp — idry, 

and write all the derivatives to the left we get the Wheeler-DeWitt equation (see Refs. 

EE) 

[-d a d p {aP + 7) - <9 a «9 7 (a7 + (3) - dpd^ + a) + iAd^d^] * = 0. (2.16) 
If we quantize using 



A,B,G^ A, B, G, 1 ' ' 



we find (See Ref. W) 



[AB{-d A d B + id G ) +AG(-d A d G + id B ) +BG{-d B d G + id A ) +AABG] # = 0. 

(2.18) 

The equivlent of this equation for supergravity has been considered in Ref. [ (Tcfl . 
Finally we consider the new Hamiltonian constraint (2.10). Using the quantizations 
introduced above we find, respectively, 

{-2 [(a{3 + i)d a d p + (cry + P)d a d 1 + (/3 7 + a)fy0 7 ] - 2iA^SgS^+ 

+ dj + d 2 ) - + + did*)} * = 0, 

(2.19) 

and 

{2 [AB(-d A d B + id G ) + AG(-d A d G + id B ) + BG(-d B d G + id A )\ + onN 
+2(A 2 + 5 2 + G 2 ) - + ^ + ^) } * = 0, (2 - 20) 

where in (2.19) we have multiplied by A 2 B 2 G 2 to avoid derivatives in the denom- 
inators. Equations (2.17) and (2.19) are ninth order partial differential equations 
(PDE's) (eighth order if we do not include the cosmological constant term); they 
are quite complicated and will have a number of "spurious" solutions introduced by 
multiplying by sixth order polynomials. Their solutions must be related to those of 
(2.14) and (2.20) by Fourier transform. Equations (2.17) and (2.18) were studied by 
Kodama||. In this paper we will concentrate on the discussion of the solutions to 
Eq. (2.20) and their relationship with solutions to the other equations. 
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In order to study equations (2.14), (2.16), (2.18), and (2.20) we will write * = 
We~ s where W and S are functions of (M, N, L), (a, ft, 7), (A, B, G), and (A, B, 
G) respectively. In this way we get the equations shown in Appendix A. We will 
look for particular solutions having the property that S satisfies a Hamilton- Jacobi 
equation (much in the spirit of the WBK approximation scheme) In this way each 
of the equations in Appendix A divides into two; an equation for S and another for 
W. The Hamilton- Jacobi equations obtained (corresponding to (2.14), (2.16), (2.18), 
and (2.20) respectively) are 



-2 [MNd M Sd N S + LMd L Sd M S + LNd L Sd N S] + 



+ 



0. 



(2.21) 



C(M 2 + N 2 + L 2 ) - ±C + + M L ¥ 1 ) - 2A MNL 

(a(3 + ^)d a Sd p S + (cry + ^Sd^S + (/3 7 + ot)d p Sd-,S + iA d a Sd p Sd^S = 0, (2.22) 

ABd A Sd B S + AGd A Sd G S + BGd B Sd G S+ , . 

+i(AB3 G S + AGd B S + BGd A S) - A ABG = 0, { ] 

ABd A Sd B S + AGd A Sd G S + BGd B Sd G S+ 

+i(ABd G S + AGd B S + BGd A S) - (A 2 + B 2 + G 2 )+ ( 2 .24) 
±(^f + ^£ + ^)+AABG = 0. 

We will now discuss some solutions (among them several new ones) to the Hamilton- 
Jacobi equations (2.21-2.24). Starting from (2.21), and putting A = 0, we will take 
an ansatz of the form 

S = a(^ + ^. + Mfy bL + c M + d N + e, (22,) 

where a, b, c, d, e are constants that we have to fix. The Moncrief-Ryan solution[0] 
to (2.21) is contained in this family and corresponds tob = c = d = e = 0. Notice 
that the addition of linear terms is not trivial due to the non-linear character of the 
Hamilton- Jacobi equation. Substituting (2.25) in (2.21) we find that the equation is 
satisfied if the constants in (2.25) are solutions of the following equations 



-4a 2 + C = 0, 
be + 2ad = 0, 
bd + 2ac = 0, 
cd + lab = 0, 



(2.26) 



and e is arbitrary. For Lorentzian signatures (( = +1) the possible solutions to (2.26) 



arc 



(fl,6,c,d) = (|,0,0,0), (|,-1,-1,-1), (|,1,1,-1), (|,1,-1,1), 

(-|,0,0,0), (-1,1,1,1), (-1,-1,-1,1), (-|,-1,1,-1), (-±,1,-1,-1). 

(2.27) 

If £ = — 1 (Euclidean signatures) the solutions to (2.26) are given by 

(a,b,c,d) = (|,0,0,0), (|, -i, -i, -i), (f, i, i, -i), (|, i, -i, i), 
(-f, 0,0,0), (-|,-i,-i,z), (-§,-i,z, 

(2.28) 
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The solutions shown above seem to be in correspondence with the analytical solutions 
known in closed form for the Bianchi IX model[[Ll|]. It is possible that there are 
no more solutions for S that can be written in analytic form[|T2|. The difference 
between the solutions for Euclidean and Lorentzian formulations is the appearance of 
a global, purely imaginary, factor. We have not been able to find solutions for non-zero 
cosmological constant. Equation (2.22) has been studied in detail by Kodama0. Here 
we give the known solutions for completeness (in this case the cosmological constant 
A must be different from zero) 



b ~ 2A 



a 2 + (3 2 + Y + 2a/3 7 . (2.29) 



We are not aware of any solution for A = 0. 

Equations (2.23) and (2.24) differ only in the potential term and the sign of the 
term with the cosmological constant. For A = we try solutions of the form 

(ad \r* RC\ 
— + — + —j+bA + cB + dG + e. (2.30) 

It is easy to verify that (2.30) are solutions to the Hamilton- Jacobi equation (2.23) 
provided that the constants satisfy the following conditions (e arbitrary) 



a 2 + ia = 0, 
be + id + lad = 0, 
bd + ic + 2ac = 0, 
cd + ib + 2ab = 0, 



(2.31) 



the solutions to the previous equations are 



(a, b, c, d) = (0,0,0,0), (0, -i, -i, -i), (0,i,i,-i),(0,i,-i,i), (0,-i,i,i), 
(-1,0,0,0), 

(2.32) 

If we consider instead equation (2.24) with its potential term, (2.30) is a solution (for 
A = 0) when the constants a, b, c, d, e are solutions to the equations 



(2.33) 



2a 2 + 2ia - 1 = 0, 
bc + id + 2ad = 0, 
bd + ic + 2ac = 0, 
cd + ib + 2ab = 0, 

where, as before, e is arbitrary. The solutions to (2.33) are 

(a,b,c,d) = (if, 0,0,0), (if, -1,1,1), (if ,1,-1,1), (if,l, 1,-1), 
(if, -1,-1,-1), (^i,0,0,0), 0=i,l,l,l), (^,1,-1,-1), (2.34) 
(if, -l,l,-l),(i=i, -1,-1,1). 

As before we know of no solution for the A ^ case. 
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Ill Discussion of Solutions 



In order to display the solutions given in the previous section we will give them in 
terms of the Misner parametrization of the Mixmaster model ||. For the Lorentzian 
ADM equations (2.1) the metric variables given in Eq. (2.2), M, N, L, are, in terms 
of the Misner variables, 



M = e 2a e 2f3 \ N = e 2a e~ l3 + + ^-, L — e 2a e ~P+-^P- _ (3.1) 

This means that the Lorentzian solutions for S given in Eq. (2.25) have the form 
(with the trivial constant e = 0) 

S = e 2a {a[e-^+ + 2e 2f3 + cosh 2y/3/3-] + ce 2/3+ + 

+e- /3+ {(6 + rf)coshv / 3 / 9_ + (d-b)smhV3{3^}). (3.2) 

Notice that the solutions (2.30) for the real Ashtekar equations have exactly the same 
form as (3.2) since A = M, B = N, G = L, but with different constants a, • • ■ ,d. 

We would like to display the solutions for a = const, (average radius of the 
universe constant) in the /3 + /3_-plane. The solution for S (3.2) with b = c = d = is, 
as was mentioned above, the solution given by Moncrief and Ryan|§ (and by Graham 
as the bosonic sector of supergravity||) and Figure 1 is a three-dimensional plot of 
|\l/| 2 = e~ 2S for a = in the (3 + (3- -plane. The line shown on the plot represents a 
contour of the potential V{(3±) given in Ref. which drives the Mixmaster model. 
The potential has a triangular symmetry, where rotation by 7r/3 in the /3 + /3_-plane 
leaves V invariant. There are also soft "channels" where V goes to zero for large 
values of that begin at the corners of the triangle and run directly to infinity. 

These channels become exponentially narrower at large values of (5 + and The 
straight lines that define the center of these channels each represent the Taub model, 
a special case of the Mixmaster model|0|. The point f3 + = (3- = is the k = +1 



Robertson- Walker universe. The solution shown in Fig. 1 is peaked over (3+ = (3- = 
and has a roughly triangular form with the points of the triangle lying in the directions 
of the three "Taub" channels of the potential. 

At first glance there seems to be a large number of new solutions with b, - ■ ■ ,d 
nonzero, but these solutions share the triangular symmetry of the Misner potential, 
so some of them are just copies of the others related by a rotation by 7r/3 in the /3+/3_- 
plane. Of the Lorentzian "wormhole" [13 solutions of Eqs. (2.26), with a = 1/2, there 



is one isolated new solution with b = c = d = — 1, which has the triangular symmetry 
of the potential. This solution, for a = is shown as a graph of e~ 2S in Fig. 2. It 
has a peak over (3+ = (3- = as before, but it seems to single out the Taub model 
channels with "arms" where the solution goes asymptotically to one as the distance 
out along the three channels becomes infinite. The "arms" become narrower rapidly 
as the distance from f3 + = (3_ = becomes large. For the arm along the /5 + -axis at 
a = and for small /?_ and large (3 + , e~ 2S becomes 

exp{-12e 2/3+ /3!}, (3.3) 
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a Gaussian in /3_ with width yl/12e~ /3+ . We will discuss the problem of "normaliza- 
tion" of these solutions below. 

Of the last three Lorentzian a = 1/2 solutions given in (2.17), only one is relevant, 
with the others given by 7r/3 rotations in the /3 + /3_-plane. We choose b = d = 1, 
c = —1, and for these constants e~ 2S is shown in Figure 3. This solution is very 
strongly peaked over almost the entire Taub model line that is represented by /3_ = 0. 
As before, for large /3+ and small f3- we have have exactly the same form for e~ s as 
given in (3.3) and e~ 2S — > 1 for (3- = and (3 + — > oo. This solution falls rapidly to 
zero for (3+ < 0, so the peak near [3 + = (3- = evident in the solutions given above 
has disappeared. 

For the solution given by (2.30), as we have mentioned, Eq. (3.2) still describes 
the "real Ashtekar" solutions with a, ■ ■ • ,d given by (2.34). Since the only difference 
between this solution and the Lorentzian solutions is that a becomes complex, but 
with real part the same as for the Lorentzian solutions, so |\l/| 2 is the same as in the 
Lorentzian case, and Figures 1-3 give this function also. 

For the ADM Euclidean solution and the true Ashtekar solutions S becomes pure 
imaginary, and |\I/| 2 is one, so we give no graphs of these functions. Of course, there is 
no reason to suppose that |\I/| 2 in the Ashtekar case has any intrinsic meaning (such 
as the "probablility" of finding the universe at some /3 + /3_ point at a = 0), since 
there is no agreement on probability measures for these complexified theories. 

It is obvious, however, that our graphing of |\I/| 2 implies that we are thinking 



of the Hartle-Hawking [jig] definition of the probablility associated with the wave 
function of the universe. The definition of probability measures on solutions to various 
formulations of quantum gravity has been a difficult problem. The Hartle-Hawking 
definition for ordinary ADM variables is one possibility, but other definitions in terms 
of superspace currents have also been proposed ||16||. For complexified variables of the 
Bargmann- Segal or Ashtekar type, the construction of probability measures is even 
more complicated, so we will not attempt to give any such measures for the solutions 
given in (2.22) and (2.23). 

The situation is slightly different for the ADM solutions (2.25) and the real 
Ashtekar-like solutions (2.30). Notice that the derivative terms in (2.14) form a 
Laplace-Beltrami operator for the superspace metric 
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NL 


L 2 



(3.4) 



so we can define a conserved superspace current 



Jo = 



/MN^*d N ty + ML^*d L ^ - MN^d N ^* - ML^d L ^* 

MN^*d M ^ + NL^*d L ty - MN^d M ^* - NL^d L ^>* \ (3.5) 
V ML^*d M ^ + NLV*d N V - MLtyd M ^* - NL^d N ^* 



valid for any solution to (2.14). Unfortunately ^ = We~ s for W = const, and S 
given by (2.25) is not a solution to (2.14). For b = c = d = Moncrief and Ryan[Q] 
showed that the factor ordering for the derivative terms of (2.14) that allows W = 
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const, is 

jd M MNLd N + ^d L MNLd M + -^d N MNLd L , (3.6) 

and a conserved superspace current for this equation is 

j k = {MNL)jI (3.7) 

where j§ is defined above. Unfortunately both these currents are zero for ^ given 
by (2.25) with b = c = d = 0, since \I/ is real. It might appear that the solutions 
(2.30) would have nonzero currents since they are complex, but the equation of motion 
(2.20) has pure imaginary terms, and we have been unable to find a conserved current 
that is compatible with it, so we cannot say that such a current would be nonzero. 
The only interesting "probability" measure is then the Hartle-Hawking |\1/| 2 which 
we have given in Figs. 1-3. 

We might worry about the "normalization" of the new functions given by (2.25) 
with b, c, d since they do not fall off rapidly for large f3± (we are thinking 
of "normalization" in the sense of an ADM equation solution with a taken as an 
internal time, which means that / d/3 + d(3_ \^\ 2 over the /5 + /3_-plane should be finite). 
However, the "arms" of these solutions that remain finite as we move out the Taub 
channels begin near /3 + = /?_ = and become very narrow rapidly, and the fact that 

roc roo 1 j 

/ dp+ d[3_exp{-12e 2/3 +p 2 _} = -Jn/3 (3.8) 

JO J-oo 2 v 

implies that J df3 + df3-\^\ 2 will remain finite. 



IV Conclusions and Suggestions for Further Re- 
search 

In this paper we have studied the use of canonical transformations to find particular 
solutions to quantum gravity, at least in the minsuperspace sector. What we showed 
was that any solution, S, to the Einstein-Hamilton- Jacobi equation can be used to 
generate a canonical transformation that leads to a solution of the form \1> = We~ s . 
The prefactor W serves to allow us to adjust the factor ordering of the Wheeler- 
DeWitt equation, since any choice of W is valid for some factor ordering. For some 
solutions S the factor ordering that allows W = const, is relatively simple, as in (3.6), 
but for others, the only such factor ordering we have found is somewhat clumsy (see 
Appendix B), but simpler expressions for W probably exist. 

We also showed that there are, in principle, a large number of functions S that 
are solutions to modified Einstein-Hamilton- Jacobi equations that represent other 
canonical transformations that do not annihilate the curvature term in the Hamil- 
tonian constraint. Each of these canonical transformations alows us to find exact 
solutions of the form We~ s . The problem is that one would like to have an analytic 
solution for S, and, in principle, one must solve a nonlinear functional differential 
equation such as (1.5) in order to find S. If one does not insist that S be a true solu- 
tion to the Einstein-Hamilton- Jacobi equation, but rather a solution to the modified 
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equation that would come from (1.9), it may still be possible to simplify the curvature 
term in the equation to the point where S may be found easily. 

Seen from this viewpoint, the Ashtekar variables, based on a complex canonical 
transformation, have the advantage that S is known exactly for the transformation, 
and one can construct solutions of the form We~ s with W = const, easily. Unfor- 
tunately there are two drawbacks to this procedure. One is that often the S one 
calculates is zero (for example, S = is a possible solution to Eqs. [2.22] and [2.23] 
with A = 0). The other is that since the transformation is complex, there is a need 
for a "reality condition" on the functions \&. Because of this we investigated a real 
transformation similar to the one that generates the Ashtekar variables that requires 
no reality condition, but makes it necessary to solve a more complicated equation for 
S. 

In order to give concrete examples of the idea of finding solutions given by gener- 
ators S obtained by solving the equivalent of the Einstein-Hamilton- Jacobi equation 
we invesigated Bianchi-type minisuperspaces, and by writing the equations in terms 
of variables of the Ashtekar type, we were able to find several solutions for S, some of 
them new. Unfortunately, the only explicit factor ordering which allowed W = const, 
we were able to find is not very appealing. 

Of course, minisuperspace solutions, while they may be of interest as "wave func- 
tions of the universe", are perhaps better thought of as clues to finding solutions to 
full quantum gravity. The Chern-Simons solution in quantum gravity was found in 
just this way. In the future we might hope to find solutions to full quantum gravity 
that are suggested by the minisuperspace exact solutions given above. This would re- 
quire solutions S to the functional Einstein-Hamilton- Jacobi equation or its analogues 
where the curvature term is not annihilated. 

This idea is perhaps the strongest suggestion for further research. It has proved 
fruitful to look for solutions to the Einstein-Hamilton- Jacobi equation and its ana- 
logues, since it has been possible to find a number of exact solutions to a problem 
such as the Mixmaster model which has a reasonably complicated structure, and one 
could easily have assumed that it would be impossible to find any analytic particular 
solutions to the problem. That it is fairly easy to find such solutions in this min- 
isuperspace case leads one to believe that it would not be impossible to find such 
solutions in full quantum gravity. The existence of the Chern-Simons solution seems 
to point in this direction. 
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A Equations for S and W 



The equations obtained by writing the wave function as ^ = We s and substituting 
it in (2.14, 2.16, 2.18, 2.20) are 

-2 [MNd M d N W + LMd L d M W + LNd L d N W- 
-MN(d M d N S)W - LM(d L d M S)W - LN(d L d N S)W- 
-MNd M Sd N W - MLd M Sd L W - LNd L Sd N W- 
-MNd M Wd N S - MLd M Wd L S - LNd L Wd N S+ 
MN{d M Sd N S)W + ML(d M Sd L S)W + LN(d L Sd N S)W] + 



C(M 2 + iV 2 + L 2 ) - \C + ^ + ^f) - 2AMNL 



W = 



-(a(3 + 7) [d a d p W - (d a d S)W - d a Sd p W - dJVdpS + (d a Sd p S)W] - 
-(cry + /3) [<9 a <9 7 W - (<9 a <9 7 S)V - <9 al S<9 7 W - cW<9 7 S + (daSd^W] - 
-(07 + a) [c^W - (SL^S^V - d^SdpW - d^WdpS + (a^fyS^W] + 
+iA [d a d p d y W - (d a dpd 7 S)W - (d a dpS)d 7 W - (d a d^S)d p W- (A.2) 
-(dpd 7 S)d a W - (d a dpW)d 7 S - (d a d 7 W)d p S - (d f3 d J S)d a W+ 

+(d a d l3 s)(d 1 s)w + {d a d 1 s){dpS)w + (d f3 d 1 s)(d a s)w+ 

+d a SdpSd 1 W + d a Sd^Sd p W + dpSd^SdaW - {d a Sd p Sd y S)W] = 

-ABd A d B W - AGd A d G W - BGd B d G W+ 
+ABd A Wd B S + AGd A Wd G S + BGd B Wd G S+ 
+ABd A Sd B W + AGd A Sd G W + BGd B Sd G W+ 

+AB(d A d B S)W + AG(d A d G S)W + BG(d B d G S)W- (A3) 
-AB(d A Sd B S)W - AG{d A Sd G S)W - BG(d B Sd G S)W+ 
+i (ABd G W + AGd B W + BGd A W) - 

-i {AB(d G S)W + AG{d B S)W + BG{d A S)W) + A ABG W = 

-ABd A d B W - AGd A d G W - BGd B d G W+ 

+ABd A Wd B S + AGd A Wd G S + BGd B Wd G S+ 

+ABd A Sd B W + AGd A Sd G W + BGd B Sd G W+ 

+AB(d A d B S)W + AG(d A d G S)W + BG{d B d G S)W- 

-AB(d A Sd B S)W - AG{d A Sd G S)W - BG(d B Sd G S)W+ ( AA > 

+i (ABd G W + AGd B W + BGd A W) - 

-i (AB(d G S)W + AG(d B S)W + BG(d A S)W) - A ABG W+ 
+{A 2 + B 2 + G 2 )w _i^ + A^ 1 + ^ W = Q 

The equations for W are then 



MNd M d N W + LMd L d M W + LNd L d N W- 
-MN(d M d N S)W - LM(d L d M S)W - LN{d L d N S)W- 
-MNd M Sd N W - MLd M Sd L W - LNd L Sd N W- 
-MNd M Wd N S - MLd M Wd L S - LNd L Wd N S = 
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-(a(3 + 7) [d a dpW - (d a d p S)W - d a Sd p W - d a Wd p S] - 

-(cry + 0) [d^W - (d a d 1 S)W - daSd^W - d a Wd 7 S} - 

-(07 + a) [d^dpW - (d^dpS)W - d^Sd p W - d 7 Wd p S} + 

+iA [d a d p d 7 W - (d a dpd^S)W - (d a d p S)d 7 W - (d a d 7 S)d p W- (A.6) 

-(d fi d y S)d a W - {d a d l} W)d 1 S - (d a d^W)d p S - {dpd y S)d a W+ 

+{d a d p s){d 1 s)w + (d^s) (dpS)w + {d p d 1 s){d a s)w+ 

+d a Sd p Sd^W + d a Sd^Sd p W + dpSd^SdaW] = 

-ABd A d B W - AGd A d G W - BGd B d G W+ 

+ABd A Wd B S + AGd A Wd G S + BGd B Wd G S+ 

+ABd A Sd B W + AGd A Sd G W + BGd B Sd G W+ (A.7) 

+AB(d A d B S)W + AG{d A d G S)W + BG{d B d G S)W 

+i (ABd G W + AGd B W + BGd A W) = 

where in each case S must be a solution for the corresponding Hamilton- Jacobi equa- 
tion. Notice that the equations for W that come from (A. 3) and (A. 4) appear to 
be the same, but the functions S that appear in them are different because they are 
solutions to different Hamilton- Jacobi equations. 



B Factor ordering for the ADM solutions 

If we want W = const., one possible factor ordering for the derivative terms in (2.14) 
is 

-2 ( MN-d MPl d N + ML-d M P2d L + NL-d N p 3 d L ) , (5.1) 

V Pi P2 P3 J 

where pi, p 2 , p 3 are functions of M,N,L. Inserting \1> = W e~ s into (2.14) with the 
factor ordering (B.l) and W = const, gives (for S given by [2.25]) 

I MN LN LM\ aMdp x / LM LN MN\ 
a \L~ + ~M + ~WJ ~ ~p^~dM \ W + W + ~L~) ~ 

ah dp 2 ( LM LN MN\ aN dp 3 /ML NL MN\ 
~~p~2~~dL \~N~ + ~L~) ~ ~~frdN \~N~ + W ~ ~L~) ~ 

pi dM p 2 dL p 3 dN y ' 

This is one equation for the three unknowns pi, p 2 , P3, so a solution is always possible. 
If, as an example, we take pi = NL, p 2 = MN, we find that p 3 becomes 

1#Wf (2(^ + bL-f) ( -bL/2a + MNL \ \ 

p 3 = MNL exp { ; v L =JU= tan" 1 . ' \ . 5.3) 

[ \JM 2 + (1 - b 2 )L 2 \y/M 2 + {l-b 2 )L 2 ) J 

This solution is not elegant, and there are probably more symmetric solutions that 
would be simpler. 
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FIGURE CAPTIONS 



Figure 1. The square of the wave function ^ = e~ s with S given by (2.25) with 
a — 1/2, b — c — d — e — 0. The solid line is an equipotential of the Mixmaster 
potential. This solution is peaked around the Robertson- Walker universe which is 
represented by (3 + = (3- = 0. 

Figure 2. |\I>| 2 = e' 2S for S given by (2.25) with a = 1/2, b = c = d = -1, e = 0. The 
solid line is an equipotential of the Mixmaster potential. While the wave function is 
still peaked around the Robertson- Walker universe, it has "arms" that single out the 
Taub models represented by the "channels" in the potential. 

Figure 3. |\1>| 2 = e~ 2S for S given by (2.25) with a = 1/2, b = d = 1, c = -1, e = 0. 
The solid line is again an equipotential of the Mixmaster potential. This solution is 
no longer peaked about the Robertson- Walker model, and lies almost entirely over 
the P + > portion of the Taub model represented by f3- = 0. 
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